Experimental Proof of a Magnetic Coulomb Phase 
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Spin ice materials are magnetic substances in which the spin directions map 
onto hydrogen positions in water ice. Recently this analogy has been ele- 
vated to an electromagnetic equivalence, indicating that the spin ice state is 
a Coulomb phase, with magnetic monopole excitations analogous to ice's mo- 
bile ionic defects. No Coulomb phase has yet been proved in a real magnetic 
material, as the key experimental signature is difficult to resolve in most sys- 
tems. Here we measure the scattering of polarised neutrons from the proto- 
typical spin ice Ho 2 Ti 2 7 . This enables us to separate different contributions 
to the magnetic correlations to clearly demonstrate the existence of an almost 
perfect Coulomb phase in this material. The temperature dependence of the 
scattering is consistent with the existence of deconfined magnetic monopoles 
connected by Dirac strings of divergent length. 
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Cooperative paramagnets are of topical interest on account of their diverse phases (1—4), 
phase transitions (5) and excitations (6-9), as well as serving as models for micromagnetic ar- 
rays (10), chemical structures (11) and colloidal films (12). Recent theoretical work (13, 14) 
has recognised that the cooperative paramagnetic state should differ fundamentally from the 
conventional paramagnetic state in the form of the spin correlations at large distance. In a 
conventional paramagnet the spin-spin correlation function decays like a screened Coulomb in- 
teraction \e~ Kr where k is the inverse correlation length, while in a large class of cooperative 
paramagnets, it is predicted to decay like a dipolar interaction ~ V r V r ' r r^7. This behaviour 
is a consequence of 'ice-rule' type constraints on the local spin configurations as was previ- 
ously noted for the analogous case of polarisation correlations in ice-rule paraelectrics (15). 
The pseudo-dipolar correlations in direct ('real') space Fourier transform into a set of so called 
'pinch point' singularities in reciprocal space. These features, which resemble bow ties, afford 
an unambiguous signature that may be characterized in a scattering experiment. However, al- 
though pseudo-dipolar correlations have been predicted for spin ices (13, 14), frustrated antifer- 
romagnets (13, 14, 16) and dimer systems (17-19), pinch points have only been clearly observed 
in the field-induced kagome ice phase of spin ice (20, 21), their absence in other experimental 
systems being notable (1, 22, 23). 

In the particular case of spin ices such as Ho 2 Ti 2 7 (2, 4), unpolarised neutron scattering 
has established that the dipolar spin ice model, in which the rare earth ions are coupled by the 
dipolar interaction and competing superexchange, gives a very accurate description of bulk and 
microscopic properties (24, 25). However, a projective equivalence between the long ranged 
dipolar interaction and a near-neighbour ferromagnetic interaction, in which the ice rules alone 
constrain the ground state spin configurations (14, 26), implies that the spin-spin correlations of 
Ho 2 Ti 2 7 will be very close to those of the full ice rule manifold of the near neighbour spin ice 
model (2). Despite the strong expectation of the existence of pinch points in the zero field spin 
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ice state (13, 14), none have yet been resolved (20, 24, 25, 27, 28). 

It is particularly important to establish the existence of pseudo-dipolar correlations in a real 
spin ice as Castelnovo, Moessner and Sondhi (8) have predicted that thermal defects in the 
spin ice state are equivalent to ionic defects in water ice, not only in a statistical sense (which 
follows from the spin ice mapping (2)), but also at the level of electrodynamics. In the spin ice 
ground states the magnetization M(r) and the magnetic H-field are divergence free, but thermal 
excitations create deconfined sources and sinks in M and H that may be described as magnetic 
'monopoles' (8, 9). The divergence free constraint on the coarse grained M(r) is a consequence 
of the pseudo-dipolar spin correlations, while the Coulombic attraction of the excitations is a 
consequence of physical dipolar interactions between spins. As dipolar interactions are not in 
doubt, the observation of dipolar correlations (pinch points) would therefore prove the existence 
of the postulated 'magnetic Coulomb phase' [] 

Neutron scattering, when appropriately energy integrated, estimates the scattering function 
S a ^(Q) in reciprocal space (here a, (3 = x, y, z), which is the required Fourier transform of 
the thermally averaged two-spin correlation function. Our polarised neutron experiments were 
configured to measure two independent components of the tensor S a ^(Q) that we label as 
spin flip (SF) and non spin flip (NSF), but only one of these two components (the SF) would 
be expected to contain pinch points (see supporting online text). The previous unpolarised 
studies (20, 24, 25, 27, 28) measured the sum of the SF and NSF scattering. 

A single crystal of Ho2Ti 2 Oy was mounted on the D7 polarization analysis spectrometer at 
the Institut Laue-Langevin, France, to allow us to map diffuse scattering over the entire (h, h, I) 
plane. Here x (y) is defined to be parallel (perpendicular) to the scattering vector Q in the 
horizontal scattering plane. With vertical incident neutron polarization, the SF and NSF cross 
sections yield information on S yy (Q,) and S ZZ (Q) respectively (see supporting online text ). 

'in this context we take a Coulomb phase to be a thermodynamic phase in which the elementary excitations are 
monopoles. 
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In Fig. [T]a-c we show the scattering in the SF and NSF channels at T = 1.7 K and the 
total, which would be observed in an unpolarized experiment and can therefore be compared 
to the observations and calculations of Refs. (24) and (28). Pinch points are clearly visible 
in the SF cross section at the Brillouin zone centres (0,0,2), (1,1,1) and (2,2,2), while the NSF 
scattering does not contain pinch points. The total scattering (SF+NSF) reveals the pinch points 
only very weakly, if at all (Fig. [TJ:), as the NSF component dominates near the zone centre. 
This is explicitly illustrated in Fig. [2J3, where cuts across the zone center show that the strong 
peak at the pinch point in the spin flip channel is only weakly visible in the total. The use 
of polarized neutrons thus extracts the pinch point scattering from the total scattering, and the 
previous difficulty in resolving the pinch point is clearly explained. 

The projective equivalence of the dipolar and near-neighbour spin ice models (14) suggests 
that above a temperature scale set by the differences between the two models, the scattering 
from Ho 2 Ti 2 7 should become equivalent to that of the near neighbour model. T = 1.7 K 
should be sufficient to test this prediction as it is close to the temperature of the peak in the 
electronic heat capacity that arises from the spin ice correlations (1.9 K) (24). Simulations 
of the near neighbour spin ice model are shown in Fig. [TJl-f . Here an effective J = 1.26 K 
is chosen to give the best fit to experiment. Note that it differs from the effective J c s = 1.8 
K derived from comparing specific heats (4), but there is no reason to expect these effective 
scales to be identical. Referring to Fig. [Tj the experimental SF scattering appears to be very 
well described by the near neighbour model, while the NSF scattering is not reproduced by 
the theory. However, we have discovered that S(Q) cxpenment /^(Q) theory is approximately the 
same function /(Q) for both channels. Thus, since the theoretical NSF scattering function is 
approximately constant, we find /(Q) ~ S '(Q)^5^ imcnt . This function may be described as 
reaching a maximum at the zone boundary and a finite minimum in the zone centre. Using the 
above estimate of /(Q), the comparison of the quantity 5(Q)^ enmcnt //(Q) with £(Q)s h ^ ory 
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is considerably more successful. Differences are less than 5% throughout most of the scattering 
map (see supporting online text and Fig. S[T|). 

In Fig. [2^ and b we show cuts through the pinch point at (0,0,2), at 1 .7 K. It has the form 
of a low sharp saddle in the intensity. In the limit of zero temperature, when the ice rules are 
perfectly obeyed, theory predicts that the pinch point will be a sharp bow tie, with a singularity 
at the centre: formally S yy (q) oc q±/{q 2 + q±)- Here q±, q\\ are components of the scattering 
vector Q measured from a reciprocal lattice point. In order to better resolve the lineshape of the 
pinch point we preformed an analogous polarised neutron experiment on the IN 12 triple axis 
spectrometer, which affords higher resolution than D7: to achieve this we used incident neutron 
wavevector ki = 1.5 A -1 (Ei = 4.73 meV) and 30' collimators before and after the sample. To 
compare with theory we use an approximation to an analytic expression given by Henley (13) 
(see also ref. (29)). In the vicinity of the (0,0,2) pinch point this becomes 



qf-2 + & 
qf-2 + ql + ql + £ 



S yy (q h , q k , qi) oc 2 , i T" f, c - 2 - (D 



Here £ icc is a correlation length for the ice rules that removes the singularity at the pinch 
point (29). In Fig. [2b we show the high resolution data from INI 2, scaled to those of D7. 
By convolving with the triple axis resolution function, the IN 12 data can be described by this 
form, with a correlation length £ icc « 182 ± 65 A, representing a correlation volume of about 
25000 spin tetrahedra. The correlation length rises to 234 ± 67 A at 1.3 K and has a tempera- 
ture variation consistent with an essential singularity ~ exp(£>/T), with B — 1.7 ± 0.1 K (see 
Figure 3c). The physical interpretation of £ ice is discussed below. 

As described above, the scattering in the NSF channel directly visualises corrections to the 
ideal spin ice behaviour. It is concentrated around Brillouin zone boundaries, as previously 
observed in unpolarised cross sections for both Ho 2 Ti 2 7 (28) and Dy 2 Ti 2 7 (27). In the latter 
case, this enhanced zone boundary scattering was well described by fine tuning further neigh- 
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bour superexchanges in the dipolar model, though it is not a direct consequence of these small 
terms (25). The NSF scattering (i.e. /(Q)) also shows a pronounced symmetric minimum at 
each Brillouin zone centre, which, with superior resolution, is revealed to be roughly as sharp 
as the maximum in the SF scattering. This indicates a partial suppression of magnetization 
fluctuations on the length scale of the ice rule correlation length, yet one that is insufficient to 
remove the pseudo-dipolar correlations up to that scale. Projective equivalence in Ho 2 Ti 2 7 
is thus proved to be robust to the expected r~ 5 corrections (14) and superexchange (25). The 
distinct structure of the NSF scattering persists to temperatures as high as 10 K, suggesting a 
simple and generic correction to the near neighbour model. 

In Fig. [3^ we illustrate the general effect of increasing temperature on the SF scattering 
pattern. Empty areas of S yy (e.g. near 1.5,1.5,0) are increasingly filled in as the temperature 
rises, with a thermal contribution that is independent of wavevector (apart from the Ho 3+ form 
factor), indicating uncorrelated point defects, or monopoles which remain strongly confined as 
dipole pairs. This is also illustrated in Fig. |2]1 where, with increasing temperature, there is a 
marked decrease in peak intensity and an increase in the background on which the peak stands. 
As shown in Fig.[3]i, this contribution can be generally fitted by the following form: 

Tm exp(-2J/T)+exp(-8J/r) 
1 ' 1 + exp(-2 J/T) + exp(-8 J/T) ' 1 ; 

where I(T) is the intensity and J = J cff = 1.8 K, the effective near neighbour exchange appro- 
priate for Ho 2 Ti 2 07 (4). The two terms are the cost of creating singly 'charged' (monopole) or 
doubly 'charged' thermal defects in the ice rules respectively. 

Experimental evidence for the additional contribution of deconfined defects is provided by 
the broadened pinch point scattering, discussed above. We interpret this length scale as the max- 
imum length of the effective Dirac strings that connect deconfined magnetic monopoles (8, 9). 
A Dirac string is a line of overturned dipoles connecting two monopoles, i.e. a ferromagnetic 
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fluctuation in spin ice. At high temperatures the proliferation of bound defects will both disrupt 
existing strings, and reduce the mean free path for diffusing monopoles, reducing the maximum 
length in the Dirac string network. As the temperature is reduced, the thermal defect population 
decreases and we see the correlation length of the Dirac string network diverges as approxi- 
mately exp(£>/T), as shown in Fig. [3];, with the observed value of B very close to the effective 
exchange J e g = 1.8 K (4). Such a temperature variation of the correlation length is the same as 
that of the one dimensional Ising ferromagnet, where the interpretation of the correlation length 
as an effective maximum string length is exact. 

The observed temperature dependent contributions to the diffuse scattering are thus well 
described by the sum of two terms, corresponding to confined and deconfined monopoles re- 
spectively. The integrated intensity shows that the confined monopoles are by far the dominant 
contribution at finite temperature, but the importance of the deconfined pairs increases as the 
temperature decreases, and at low temperatures they are deconfined over macroscopic (at least 
micron) length scales. This behaviour is fully consistent with a picture of a weakly dissociated 
electrolyte where the energy cost of dissociation (including screening) is small compared to the 
energy of ion pair formation. In this regard, the spin ice thermal defects are fully equivalent to 
ionic defects in water ice (but note that there is no equivalent of the 'Bjerrum defect' (30)). 

In conclusion, we have studied the spin ice Ho2Ti 2 07 using polarized neutron scattering. 
We have established the validity of projective equivalence and have quantified the corrections 
to it. These corrections are summarised by the function /(Q) that dips sharply at the zone 
centre, raising the question of whether /(Q) can be captured in a generic improvement to near 
neighbour spin ice that still obviates a full treatment of the microscopic Hamiltonian. More 
importantly, we have established that Ho2Ti 2 07 exhibits an almost ideal magnetic Coulomb 
phase, the quasiparticle vacuum for Castelnovo et a/.'s magnetic monopoles (8, 9). We have 
estimated the length of the longest Dirac strings to rise to macroscopic scales as the temperature 
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passes below 1 K. 
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Figure 1: Diffuse scattering maps from spin ice, Ho2Ti 2 07. Experiment (top row, a-c) versus 
theory (bottom row. d-f) a: Experimental SF scattering at T = 1.7 K with pinch points at 
(0,0,2), (1,1,1) and (2,2,2) etc. b: the NSF scattering, c: the sum, as would be observed in an 
unpolarized experiment (eg. Refs (24, 28)). Sharp spots are nuclear Bragg peaks and the data 
have been corrected for absorption, d: The SF scattering obtained from Monte Carlo simulations 
of the near neighbour model with 12x12x12 supercell, scaled to match the experimental data, 
e: The NSF scattering from this model is Q-independent, modulated only by the Ho 3+ magnetic 
form factor, f: The total scattering of the near neighbour spin ice model. 



11 




2 
0,0,1 



a:1.7K 



o SF 
♦ NSF 



3 

CD 




g 
CD 



^ 4 

c 

CD 

£ 2 





c:1.7K 




□ 


Total 


♦ 


NSF 





SF 



♦ NSF 

o SF 



-50 K 

20 K 

10 K 

5 K 

3.75 K 

2.5 K 
-1.7 K 



-0.6 -0.3 0.3 0.6 
h,h,2 

Figure 2: Line shape of the pinch point. Top (a): radial scan on D7 through the pinch point 
at (0,0,2) (er' is the neutron scattering cross section: see supporting online text for its precise 
definition). Second (b): the corresponding transverse scan, where the lines are Lorentzian fits. 
Third (c): higher resolution data from IN12, where the line is a resolution corrected fit to the 
pinch point form Eqn. ([T]): the resolution width of IN12 is indicated as the central Gaussian. 
Bottom (d): SF scattering at increasing temperatures (the lines are Lorentzians on a background 
proportional to the Ho 3+ form factor). 
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Figure 3: Temperature dependence of defect scattering, a: the SF scattering at 5 K. b: the 
spin flip scattering at 10 K. The rising temperature washes out the sharp pattern of 1.7 K. In 
particular, simulations suggest that in the limit of perfect ice rule enforcement, the area around 
(2,2,0) contains no magnetic scattering. As the temperature rises, this area is filled in (compare 
also Fig.[T^). c: the correlation length extracted from the high resolution data is consistent with 
an exponential divergence of Dirac string length at low temperature, d: temperature dependence 
of the scattering at (0.5,0.5,2) (the pinch point background in Fig.|2]l) and (1.5,1.5,0) (indicated 
by corresponding symbols in a). These parts of the pattern are sensitive to the bound defect 
population. The lines are fits to Eqn. [2] (there is a fitted scale factor), which describes the 
formation of single and doubly 'charged' ice rule defects (see text). 
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Materials and Methods 

Crystal growth: polycrystalline Ho 2 Ti 2 07 (space group Fd3m) was synthesized from high 
purity Ho 2 3 and Ti0 2 powders by reacting in air at 1250 °C for 3 days with intermediate 
grindings. A crystal boule was grown using a four mirror optical floating zone furnace with a 
growth speed of l-2mm/h and a rotation speed of 25 rpm under flowing oxygen gas. The single 
crystal used was cut from this boule and its quality checked by x-ray diffraction. The sample is 
approximately cylindrical with radius ks 3 mm and length ps 45 mm. 

Neutron Scattering experiments on D7: The crystal was oriented with the [110] axis vertical 
(with the boule tilted at some 20° away from vertical) and clamped to a copper holder. In order 
to minimize the incoherent scattering background, the use of glue was avoided. The temperature 
was controlled using a standard Orange cryostat. Details of the polarization analysis system of 
D7 can be found in Ref. S 1 . 

With z = [110] vertical, the (h, h, I) wavevectors perpendicular to this formed a conven- 
tional xy scattering plane. The incident neutron energy was 3.5 meV and the incident neutron 
polarization P was vertical. As in previous studies (20, 24, 25, 27, 28), the static approxima- 
tion (i.e. an appropriate integration over energy) was ideally obeyed, so the scattering function 
measured the arrangement of spins at an instant in time, averaged over an ensemble of identical 
systems. Defining x as parallel to Q and y orthogonal to x and z, the SF and NSF cross sections 
are proportional to S yy (Q) and S ZZ (Q) respectively. The tetrahedral basis of the pyrochlore 
stucture has four spins: for a general scattering vector all four spins contribute to S yv (Q), but 
only two contribute to S ZZ (Q) as the Ising anisotropy means that the other two have no compo- 
nent on z. This has the strong consequence that pinch point singularities are expected in S vy (Q) 
but not in S ZZ (Q), as only S yy will contain the full set of ice rule correlations. The cross sec- 
tion for unpolarised neutron scattering is proportional to S yy (Q) + S ZZ (Q). Other contributions 
to the scattering (nuclear Bragg, nuclear thermal diffuse, nuclear spin and isotopic incoherent 
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scattering) have a negligible effect on the results presented here, so will not be discussed further. 

We measured the non-spin flip and spin flip cross sections with the polarization parallel to z 
by rotating the crystal (fl). The data in Fig. [T] of the paper were obtained by 0.5° step rotations 
and with two positions (6) of the multi-detector banks to eliminate gaps in the scattering map. At 
higher temperatures, only one detector bank position was measured, resulting in the gap that can 
be seen in Fig. [3^-b. The statistics around the (0,0,2) pinch point were improved by a factor of 
2 by extended counting at all temperatures. The classical || — _L polarization analysis technique 
can be applied at a single value of \Q\ enabling a complete separation of scattering contributions 
in anisotropic materials, not generally possible in the multidetector xyz measurement. We 
measured this cross section at a single point at 20 K to confirm our normalization procedure. 

The data were corrected as described in Ref. SI. In brief, the instrumental background 
was subtracted, an amorphous silica standard was used to normalize the supermirror analyzer 
efficiencies and a known quantity of vanadium to normalize detector efficiencies and provide 
an absolute scale. All scattering data from D7 shown in the paper is therefore a measurement of 
the absolute partial differential cross section a' = do /dVt in barn steradian -1 (formula unit) -1 . 
Additionally the data were corrected for absorption. Path dependent attenuation coefficients 
were calculated for every point in the Vt — map using the U 5-matrix approach outlined in Refs. 
S2 and S3. After this procedure the || — _L measurement confirmed that the isotopic incoherent 
scattering was as expected for Ho 2 Ti 2 7 . The spin incoherent scattering was somewhat larger 
than calculated, but only a very small amount of hydrogenous material would be required to 
produce this effect. In Fig. la sharp intensity is visible at the centre of the pinch points at 
(1,1,1), (2,2,2) and (0,0,4). This is due to finite polarization analyzer efficiency and is noticable 
only at the position of a very intense nuclear Bragg peak. As the temperature is lowered the 
sample slightly depolarizes the beam, but we confirmed that the total intensity of the nuclear 
peaks summed over both channels is constant, and corrected for the depolarization effect by 
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comparing the flipping ratio of the Bragg peak with the ideal value for the quartz standard at the 
same position. 

Neutron scattering experiments on IN12: The crystal was mounted as described above, but 
in a dilution refrigerator insert (we do not discuss the low temperature data as depolarization 
effects become more severe). IN 12 was configured with a focussing graphite monochromater, 
supermirror bender, Heusler analyzer and open-30-30-open collimators. The convolution of 
equation 1 was achieved by using a raytracing Monte Carlo integration algorithm described in 
Ref S4. 

Monte Carlo simulations: We simulated the near neighbour spin ice model using a Metropo- 
lis Monte Carlo algorithm combining single spin and loop dynamics. The single spin flip dy- 
namics establish equilibrium at high temperature and the loop moves maintain ergodicity at 
low temperature, as described in Ref. S5 and references therein. For comparison with neutron 
scattering we averaged two simulations of 12 x 12 x 12 pyrochlore unit cells, each contributing 
200 spin configurations. We estimated the correlation time and ensured that spin configurations 
collected from each simulation were separated by at least this many Monte Carlo sweeps. The 
calculated polarized cross sections are multiplied by the form factor of Ho 3+ (Ref. S6). 

To compare directly with the experimental data, we subtract the estimated spin incoherent 
cross section from the data and scale the simulated result. When invoking f(Q), we divide 
the experimental spin flip scattering by the non-spin flip and then scale the simulated data. 
In this case it is not necessary to include the magnetic form factor as it is divided out of the 
experimental data. We found no improvement in this process by modifying the scale or baseline 
of the non-spin flip data and this is unsurprising since Fig [2^ of the paper shows them to be 
essentially identical in certain directions, except at the zone centre. The comparison requires 
the exclusion of nuclear Bragg peaks and the zone centres, where the departure of the non-spin 
flip from the spin flip creates an artificial peak in the divided data. We quantify the goodness 
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of the fit using a difference map, showing here the results for the best fitting simulation. It is 
clear that the match is good throughout the scattering map, except for a specific collection of 
wavevectors, generally transverse to a pinch point, where the approximation of a single /(Q) 
must be poor. 
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S 1: Comparison of the experimental data and Monte Carlo simulations. In the top panels the 
comparison includes f(Q), i.e the experimental spin flip scattering is divided by the experi- 
mental non-spin flip scattering, and the simulation scaled, a and b show cuts along (0, 0, I) 
and (h, h, 0) respectively. In c we show the percentage difference between the calculated and 
experimental scattering. In the bottom panels we show the direct comparison (as in Fig. 1 
of the paper) in which the simulation is simply scaled to the experimental data. Because the 
experimental data is modulated along (0,0,/), but the simulation is not, significant mismatch 
appears. 
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